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Abstract 

We consider a integro-differential nonlinear model that describes 
the evolution of a population structured by a quantitative trait. The 
interactions between traits occur from competition for resources whose 
concentrations depend on the current state of the population. Follow- 
ing the formalism of [15] . we study a concentration phenomenon aris- 
ing in the limit of strong selection and small mutations. We prove that 
the population density converges to a sum of Dirac masses character- 
ized by the solution (/? of a Hamilton-Jacobi equation which depends 
on resource concentrations that we fully characterize in terms of the 
function Lp. 
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1 Introduction 



We are interested in the dynamics of a population subject to mutation and 
selection driven by competition for resources. Each individual in the popu- 
lation is characterized by a quantitative phenotypic trait x G M (for example 
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the size of individuals, their age at maturity, or their rate of intake of nutri- 
ents). 

We study the following equation 

dMt,x) = ^ ^^Jf(t)r/,(a;)-lj u,{t,x) + M,{ue){t,x), (1.1) 

where is the mutation kernel 

M,(/)(x) = - [ K{z) (fix + ez) - f{x)) dz, (1.2) 

for a G C^(M) such that fj^zK{z) dz = 0. Among many other ecolog- 
ical situations [H] , this model is relevant for the evolution of bacteria in a 
chemostat [131 [15]. With this interpretation, u'^(t,x) represents the concen- 
tration of bacteria with trait x at time t, the function rji{x) represents the 
growth rate of the population of trait x due to the consumption of a resource 
whose concentration is and the term —1 corresponds to the decrease of the 
bacteria concentration due to the constant flow out of the chemostat. This 
model extends the one proposed in [TB] to an arbitrary number of resources. 
This equation has to be coupled with equations for the resources /j, namely 

^'^^^ 1 + f^T]i{x) Ue{x) dx' ^^'^^ 

This corresponds to an assumption of fast resources dynamics with respect 
to the evolutionary dynamics. The resources concentrations are assumed to 
be at a (quasi-)equilibrium at each time t, which depends on the current 
concentrations u^. 

The limit e — )• corresponds to a simultaneous scaling of fast selection and 
small mutations. It was already considered in [15]. The following argument 
explains what limit behaviour for can be expected when £ — >■ 0. Defining 

ipe as 

Me = e'^"/^ or (fe = e \ogUe, (1.4) 

one gets the equation 

k 

d,ip, = - 1 + H,{^,), (1.5) 

i=l 



2 



where 

H,{f) = [ K{z) (e(/(-+- - 1) dz. (1.6) 
Jr 

At the hmit e — t- the Hamiltonian simply becomes 

H{p) = [ K{z) {eP' - 1) dz. (1.7) 
Jr 

So one expects Eq. (11. 5p to lead to 

k 

dt^ = J2 W mix) - 1 + H{d,y^), (1.8) 

i=l 

for some /j which are unfortunately unknown since one cannot pass to the 
limit directly in (II. 3p . Therefore one needs to find a relation between the cp 
and the Jj at the limit. Under quite general assumptions on the parameters 
(see Lemma [3A] below), the total population mass is uniformly bounded 
over time. This suggests that max^-gK y9(t, x) = should hold true for all 
t > 0. Together with (11. 8p . this gives a candidate for the limit dynamics 
as a solution to a Hamilton- Jacobi equation with Hamiltonian H and with 
unknown Lagrange multipliers Jj, subject to a maximum constraint. The 
limit population is then composed at time t of Dirac masses located at the 
maxima of (f{t,-). 

This heuristics was justified in [TJ] in the case of a single resource (and when 
the resources evolve on the same time scale as the population), and the case 
of two resources was only partly solved. The mathematical study of the 
convergence to the Hamilton- Jacobi equation with maximum constraint and 
the study of the Hamilton- Jacobi equation itself have only be done in very 
specific cases [13 El El ES]- In fact the main problem in this proposed model 
is that the number of unknowns (the resources) may easily be larger than the 
dimension of the constraint (formally equal to the number of points where 
¥^ = 0). 

Our goal in this paper is to prove the convergence of ip^ to a solution of (11. 8p . 
where we give a full characterization of the functions Jj. Those are no more 
considered as Lagrange multipliers for a set of constraints but are given by 
the solution ip itself. The new resulting model describes the evolution of a 
population as Dirac masses and is formally well posed. 
The general problem of characterizing evolutionary dynamics as sums of 
Dirac masses under biologically relevant parameter scalings is a key tool 
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in adaptive dynamics — a branch of biology studying the interplay between 
ecology and evolution [TSl 1221 [121 E]- The phenomenon of evolutionary 
branching, where evolution drives an (essentially) monotype population to 
subdivide into two (or more) distinct coexisting subpopulations, is particu- 
larly relevant in this framework [22| [T7t [18] . When the population state can 
be approximated by Dirac masses, this simply corresponds to the transition 
from a population composed of a single Dirac mass to a population composed 
of two Dirac masses. 

Several mathematical approaches have been explored to study this phe- 
nomenon. One approach consists in studying the stationary behaviour of an 
evolutionary model involving a scaling parameter for mutations, and then let- 
ting this parameter converge to 0. The stationary state has been proved to be 
composed of one or several Dirac masses for various models (for deterministic 
PDE models, see |H |5l [HI |20l [16] , for Fokker-Planck PDEs corresponding to 
stochastic population genetics models, see for stochastic models, see [26] . 
for game-theoretic models, see [lO])- Closely related to these works are the 
notions of ESS (evolutionarily stable strategies) and CSS (convergence sta- 
ble strategies) [221 [13], which allow one in some cases to characterize stable 
stationary states gl [IH [2Q1 [lO] . 

The other main approach consists in studying a simultaneous scaling of mu- 
tation and selection, in order to obtain a limit dynamics where transitions 
from a single Dirac mass to two Dirac masses could occur. Here again, de- 
terministic and stochastic approaches have been explored. The deterministic 
approach consists in applying the scaling of (11. ip . The first formal results 
have been obtained in [15] . This was followed by several works on other mod- 
els and on the corresponding Hamilton- Jacobi PDE [5l |25]. For models of 
the type we consider here, rigorous results (especially for the well posedness 
of the Hamilton- Jacobi eq. at the limit) mainly only exist in the case with 
just one resource, see [2] and [1] (one resource but multidimensional traits). 
The stochastic approach is based on individual-based models, which are re- 
lated to evolutionary PDE models as those in [HI [20] through a scaling of 
large population [S]. Using a simultaneous scaling of large population and 
rare mutations, a stochatic limit process was obtained in [B] in the case of a 
monotype population (i.e. when the limit process can only be composed of a 
single Dirac mass), and in |9j when the limit population can be composed of 
finitely many Dirac masses. 

Finally note that the total population of individuals is typically very high, 
for bacteria for example. This is why even stochastic models will usually take 
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some limit with infinite populations. Of course, this has some drawbacks. 
In particular the population of individuals around a precise trait may turn 
out to be low (even though the total population is large). As in the scaling 
under consideration, one has growth or decay of order exp(C/£), this is in 
fact quite common. One of the most important open problem would be to 
derive models that are both able of dealing with very large populations and 
still treat correctly the small subpopulations (keeping the stochastic effects 
or at least truncating the population with less than 1 individual). 
There are already some attempts in this direction, mainly proposing models 
with truncation, see pi] and very recently [23] . For the moment however the 
truncation is of the same order as the maximal subpopulation. 

In order to state our main result, we need some regularity and decay assump- 
tions on the rji, namely 

k 

Vi > 0, 3r]e Co(R), Vx, + + l^^'(^)l) ^ ^(^)' (1-9) 

i=l 

where Co(M) is the set of continuous function, tending to as x — )■ ±oo. 
In order to characterize the resources Ii{t) involved in (11. 8p . we introduce a 
sort of metastable equilibrium. To this aim, we need an assumption on the 
number of possible roots of the reproduction rate 

k 

3k < k, V/i . . . /fc G [0, 1], the function Jj rii{x) — 1 has at most k roots. 

i=l 

(1.10) 

We also require an invertibility condition on the matrix i]i{xj) 

Vxi . . .Xfc distinct, the k vectors of coordinates {rii{xj))^^i i are free 

(1.11) 

Then we may uniquely define the metastable measure associated with a set 
u by 

Proposition 1.1 For any dosed u C M, there exists a unique finite nonneg- 
ative measure /i(w) satisfying 
i) suppfi C u 

a) denoting = 1/(1 + / dfi{x)), 

k k 

Ii{fi) rji{x) — 1 < m cj, Vii^) — 1 = on suppfx. 

i=l 1=1 
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Now the limiting Jj are directly obtained by 

m = iMMt,-) = o})). (1.12) 

We prove 

Theorem 1.1 Assume K e C^{R), J^zK{z)dz = 0, ([L9]), ffLTOj) . ffTTTj) . 

that the initial data ^^(t = 0) > or (fs{t = 0) are C"^ , satisfy 

sup / Us{t = 0, x) dx < OO, sup ||9a;(y9e(t = 0, ■)||loo(]r) < OO, (1-13) 
£ Jr £ 

inf inf dxxfit = 0,x) > — oo, (1-14) 

£ XgR 

anc? i/iai ^^^(t = 0, ■) converges to a function ip^ for the norm || • ||y^i,oo^]]^^. 
Then up to an extraction in e, ips converges to some cp uniformly on any 
compact subset of [0, T] xM and in W^'''^{[0, T] x K) for any T > 0, p < oo and 
any compact K . In particular, (p is continuous. The function If converges to 
li in L^([0, T]) for any T > 0, p < oo, where Ii is defined from (p as in U.l^) . 
and li is approximately right- continuous for all t > 0. The function (p is a 
solution to (11. 8p almost everywhere in t, x with initial condition (p{t = 0, ■) = 
Moreover if one defines ijj = ip> — ^*L]^ /J /^(s) (is ?7i(x), then ip is a 
viscosity solution to 



dtip(t, x) 



H (^dx^ + Y,j^I,{s)ds^,{x)^ . (1.15) 



We recall that a function / on [0, +00) is approxiamtely right-continuous at 
i > if t is a point of Lebesgue right-continuity of /, i.e. 



hm- r^\f{e)-fit)\do = o. 

■'^0 s Jt 



Notice that, under the assumptions of Theorem 11.11 (psit = 0,a;) — >■ —00 
when X — )■ ±cx3 since /jj < 00 and ip^ is uniformly Lipschitz. Be also 
careful that we assume \\(peit = 0, ■) — vi/i.°°(m) — ^ even though ^^^(t = 0) 
(and thus y?") is not bounded. 

In the proofs below, C denotes a numerical constant which may change from 
line to line. 
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2 Proof of Prop. 13 



2.1 Uniqueness 

Assume that two measures /^i and /i2 satisfy both points of Prop. Il.li We 
first prove that they induce the same ressources Jj and then conclude that 
they are equal 

1st step: Uniqueness of the li. The argument here is essentially an adaptation 
of EU]. First note that 



Ji(/i2)r/,(x)-l rf/ii <0, (2.1) 



Jr V^=i / \i=l 



since /xi and fi2 are non negative and by the point ii, Yl'^=i ^iif^j) — 1 is 
non positive on u for j = 1, 2. 

On the other hand since ^JL^ ^ii^^j) vanishes on the support of fij, one 
has for instance 

k 

= 5^(/.(/il)-/.(/i2)) (l/J.M - 1), 



by the definition of /j(/U2). 
Since one has 



one deduces from (12. ip that 

Um)=UlX2), i = l...k. (2.2) 

2nd step: Uniqueness of fi. It is not possible to deduce that fii = fi2 directly 
from (12.21) . This degeneracy (the possibility of having several equilibrium 
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measures, all corresponding to the same environment) is the reason why we 
require additional assumptions on the rji. 

First of all by Assumption fll.lOp . point i and thanks to fl2.2p . we know that 
fj,i and /i2 are supported on at most k points {xi, . . . ,xi}, which are the 
roots of J2i ^iif^i)Vi{^) ~ 1 = J2i ~ 1- Therefore one may write 

k 
1=1 

Now (12. 2p tells that / rjidfii = J riidfi2 which means that 

k k 

"^aj r]i{xi) = ^afr]i{xi), \/i = l...k. 
1=1 1=1 

To conclude it remains to use condition fll.lll) and get that a} = af. 
2.1.1 Existence 

The basic idea to get existence is quite simple: Solve the equatioE0 

dtiy= ^Jiiiy)r],{x)-1^ z/, (2.3) 

and obtain the equilibrium measure fi as the limit of i/(t) as t — +00. 
This is done by considering the entropy 



/I /* /* 

L(z/) = J]log/i(z/) + dp = -^\og{l+ r],du^ + d 
i=i J i=i 



V. (2.4) 



As — log is convex and r/, > 0, then L itself is a convex function of v. 
Moreover if vii) solves (12. 3p . one has 

jL{v{t)) = -J ^^7,(z.)77.(x) - du. (2.5) 



^Existence and uniqueness are trivial for ()2.3p . for example by Cauchy-Lipschitz theo- 
rem in the set of finite positive measures equipped with the total variation norm. 



8 



Therefore one expects the hmit of v and the equihbrium measure we are 
looking for to be the minimum of L. 
Since the r]i are bounded, one finds 

L{v\ > -C + c \ dv. 



for two numerical constants C and c. Consequently L is bounded from below 
on M\{l>S) the set of nonnegative Radon measures on uj. It also gets its 
infimum on a bounded part of M\{{jS). As any ball of M\{{jj) is compact for 
the weak-* topology (dual of continuous functions with compact support), L 
attains its infimum, or 

Mo = {z/ G M\{uj), L{y) < L{u') \/u G M^(cj)} ^ 0. 

Now take any /x G Mq then take z/ the solution to (12. 3p with uit = 0) = /i. 
L(z/) is non increasing and since it is already at a minimum initially, it is 
necessarily constant. By (12.51) . this means that 

k 

liifJ') ?7i - 1 = on supp /ij. 

i=l 

Hence /j, is in fact a stationary solution to (12.31) and it satisfies point i and 
the second part of point ii of Prop. 11.11 Note by the way that the uniqueness 
argument in fact tells that there is a unique element in Mq. 
It only remains to check the first part of point ii. By contradiction assume 
that there exists a point xq G s.t. 

k 

liin) r]i{xo) - 1 > 0. 

i=l 

Let a > and define i/q = yU + aS^o G M^{u). Now compute 
L{i^a) = J dfi + a -'^log ^1 + J r]idfi + ar]i{xo)^ 




i=l 

k 



i=l 



k \ 

J2Wv^ixo)-l]+Oia''). 
i=l J 
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Thus Li^Ua) < L{fi) for a small enough which is impossible as /i is an absolute 
minimum of L. 

Consequently the first part of ii is satisfied and the proof of Prop. 11.11 com- 
plete. 



3 Proof of Theorem 1.1 



3.1 A priori estimates for Eq. ( 11.51 ) 



We denote by BViod^) ^^e set of functions on M with bounded variation on 
any compact subset of M, by M^{u) the set of signed Radon measures on the 
subset a; of M equipped with the total variation norm. 
We show the following estimates on the solution to (11.51) 

Lemma 3.1 Let (fi; be a solution to (11.51) with initial data (p^ such that 
Ir e'^^^^^/^ dx < oo, dx^^ G L°°(R) and dxx^l uniformly lower hounded. Then 
for any T > 

||c^tV^e||L°°([0,T]xR) + ||f^xV^e||L°°([0,T],Byioc(ffi)nL°°(R)) + ||f^teV^e||L°°([0,T],Mi) < C't, 
Vt < T, X G M, dxxVe{t, x) > -Ct, He{^,) > -CtS, 



'^t<T, I Ueit^x) dx < Ct, ip eit,x) < Ct£ \ogl/e. 
Jr 

where Ct only depends on the time T, J Us{t = 0)dx, \\dx(p%L°°{R) (ind 
infxdxx'-p'^ei^)- 

Proof. We start with the easy bound on the total mass. 



Step 0: Bound on the total mass. First notice that because of (11.91) . there 
exists > s.t. 

k 

V|x|>i?, ^rii{x)< 1/2. 

i=l 

Let ifj he & regular test function with support in \x\ > R, taking values in 
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[0, 1] and equal to 1 on \x\ > R+1. Using the fact that If{t) < 1, we compute 

/ l/j{x)Ue{t,x)dx <- [ 'll){x)Ue{t,x)dx 



H — / K{z){ip{x — ez) — ip{x))u^{t,x) dzdx 



e , 

< — — I il!{x)us{t,x) dx + C lus{t,x)dx. 
On the other hand as each rji > 0, one has for some constant C 



1 + J TjiU^dx 1 + J(l — ip)ue dx 
Therefore with the same kind of estimate 



d_ 
~dt 



'1 



— 1p{x))Ue{t,x) dx <C / Ue{t,x)dx 

Jr 

/ {1 — t/j{x)) U£{t,x) dx. 
Jr 



1 C 

+ - FT ^ 1 

£ \ 1 + J (1 — ^jMe dx 

Summing the two 

1 / n ^ ^ , 

1^1 — iIj{x)) U£{t, x) dx 



d f , , , If C 

— / uJt, x) dx <- TT-. 

dtJu ' ^ -e Vl + /(1- 



'il^)us dx 



f ip{x) Us{t,x) dx + C j Us{t,x) dx. 

2£ Jm. Jr 

Since the sum of the first two terms of the r.h.s. is negative if J is larger 
than a constant independent of e, this shows that J ^^(i, x) dx remains uni- 
formly bounded on any finite time interval. 

Step 1: Bound on d^^e- This is a classical bound for solutions to some 
Hamilton- Jacobi equations. Here we still have to check that it remains true 
uniformly at the e level. Compute 



i=l 



+ / K{z) d^Vejt, x + ez)- da:(p£{t, x) 



11 



We first observe that, as I- G [0, 1] and J2i Wi{^)\ ^ vi^) 

C 

~ e 

since K has compact support. This entails 

^ Jo 

from which easily follows that d^^^ G L°°([0, t^], M) for some > 0, which 
may (for the moment) depend on e. 

Now we use the classical maximum principle. Fix t G [0, T] such that Ce,t '■= 
\\dx<-Pe{t,-)\\L°°{R) < oo. For any X G M such that 9a;V2e(t, > supyda:ipe{t,y)- 
a, where the constant a > will be specified later, we have 

Jr £ \ e / 

Therefore, choosing a = ee"*" we obtain 

^ snpd:r(Pe{t,x) < C, 

for a constant C independent of e. Using a similar argument for the mini- 
mum, we deduce that > T and that dx^Pe is bounded on [0,T] x M by a 
constant depending only on T and ||9i.V5£||L°°(K)- 

Step 2: First hound on Hi,{ipf,) and hounds on dt^Pe and yj^. Simply note that 
-1 < H,{y^e{t)){x) = / K{z)e ^ dz- / K{z)dz 



< j K{z) e 

Consequently, directly from Eq. (11. 5p . 

\dt'Pe\ < vix) + C, 
hence concluding the full Lipschitz bound on ip^. 



^\\\Ox'PE\\L°^(lO,T],tl) dz <■ C 
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To get the upper bound on y^^, simply note that because of the uniform 
Lipschitz bound on ^p^ 

Ve{t, y) > (pe{t, x) -CT\y-x\, 



{t, y)dy> / e^^(*'^')/" e"^^ l^"^!/" dyyCrS e^^(*'^)/^ 



so 



Hence the bound on the total mass yields that ips < Ct£ logl/e. 

Step 3: BV hound on dx^e- As for dx^ei we begin with a maximum (actually, 
minimum) principle. First from (11. 5p 

-/ A , /" r^/ A dxxVei't.X + ez) - dxx<-Pe{t, x) 

OtOxx^e > -r]{x) + / K{z) e ^ dz 

Jr ^ 

f . ^e{t,^+ez)-^,(t,^) (d^^Jt.x + ez) - d^(pjt.x))'^ 

+ K[z)e s dz. 

Jr ^ 

The last term is of course non negative and so with the same argument as 
before, we get 

— Mdxx^e{t,x) > -C, 

dt X 

where C does not depend on e. This proves the uniform lower bound on 
dxx^e- On the other hand, for any measurable subset A of [xi,a;2], 

{IxeA - ^x<^A)dxxV^e{t, X)dx = / dxx<fe{t, x) dx 

Xl J X\ 



-2 / dxx^e{^,^)\<i,Adx 



x\ 



< dx(Peit,X2) - dx^Psit, Xl) + C \X2 - Xi\ < 2\\dx(Pe\\L^{[0,T],R) + ^ \X2 - Xi\. 

This indeed shows that dxx^eit, ") belongs to M^{[xi,X2]) with total variation 
norm less than 2||(9^y9£||2,oo([o,r],R) + C \x2 — Xi\. Thus, dxx'^e belongs to the 
space L°°([0, T],M\R)), which entails dx^fe e L°°([0, T], 5yioc(K)). 
Finally, differentiating (II. 5p once in x, one has 

\0,x^^it,x)\ < ,ix) + / ir(z)el-Nia...ll.^ l^-^-(^-^ + ^-J-^-^-(^--)l ,, 



<f]{x) + C I K{z) I \dxx^e{t,x + ee)\dedz. 
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Integrating, by Fubini 

/ \dtx^e\dx< / f]{x)dx + C / \d^^ipe\dx, 

J xi J xi J xi—ep 

where p is such that the support of K is included in the ball centered at of 
radius p. This ends the proof of all the bounds on the derivatives of ipe- 

Conclusion. It only remains to show the sharp lower bound on He{ipe)- Let 
us write 

H^{'Pe)> [ K{z)exp( [ zd^pe{t,x + eze)de] dz - ! K{z)dz. 

The BV bound on dx^Pe shows that this function admits right and left limits 
at all X G M. Let us denote dxPeit, x"*") the limit on the right and dxPeit, x~) 
the limit on the left. As d^xPe is bounded from below, we know in addition 
that 

Vx, dxPe{t,X^) >dxPei't,X~). 

By differentiating once more 
1 

zdxPe{ii X + 6ze) dO > zd^Peit, x^) 



+ I z ! ezedxxPe{t,x + e'eze)de'de 
Jo Jo 

> zdxPeit.x^) - Cez^, 
again as dxxPe is bounded from below. Finally 

Hs{Pe)> I K{z)exp{zdxPs{t,x+)-Cez^)dz- / K{z)dz 
Jr Jr 

>H{dxPe{t,x+))-Ce, 

where H is defined as in (11.71) and since K is compactly supported. Because 
we assumed that Jj^zK{z) dz = 0, we have H{p) > for any p, which ends 
the proof of Lemma 13. 1[ 
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3.2 Passing to the limit in 



From the assumptions in Theorem ll.H Lemma 13.11 gives uniform bounds on 

Therefore up to an extraction in e (still denoted with e), there exists a func- 
tion if on [0, T] X M such that dtif G L°°([0, T] x M), d^ip G L°°([0, T],BVxoc n 
L°°(M)), dtx^ G L°°([0,T],Mi^„^(M)) and d^^^ uniformly lower bounded on 
[0,T] X M, satisfying 

— > uniformly in C{K) for any compact K of [0,T] x M, 
(9j.v9e — 9j.v5 in any -^[^^([0' T],R), p < oo. 

The first convergence follows from Arzela-Ascoli theorem. For the second 

convergence, observe that \\d^^Pe\\L°°iio,T],BVi^,{M.)) + \\dtx'fe\\L°-{[o,T],M'm) ^ 
implies that d^p^ is uniformly bounded in L^{[0, T] x M) n5l^oc([0, T] x M). 
The convergence in Lf^^ follows by compact embedding. We also have < 
since otherwise the uniform bound on J^Us{t,x) dx would be contradicted. 
As the If are bounded, it is possible to extract weak-* converging subse- 
quences (still denoted with e) to some Ii{t). 
Now, we write again 

H,{^,) = J K{z) (exp(^^ zd^^,{t,x + eze)d0^ -1^ dz. 

From the L°° bound on dxfe and its strong convergence, one deduces that 

H^ip,) H{d,p) in LL- (3.2) 

Therefore one may pass to the limit in (11. 5p and obtain (II. 8p (for the moment 
in the sense of distribution; the equality a. e. will follow from the convergence 
of If in LP([0,T]), proved below). 

In addition by following [15] or [2], one may easily show that ^(t,x) = 
'p{t,x) — Yl\=i lo ^ii^) dsr]i{x) is a viscosity solution to (I1.15p . We refer the 
reader to [15] or [2] for this technical part. 

It remains to obtain (I1.12p . the approximate right- continuity of Jj for all time 
t and the convergence of If to Jj in Lp([0, T]) for p < oo. This requires some 
sort of uniform continuity on the If which is the object of the rest of the 
proof. 
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3.3 Continuity in time for the If 

First of all note that, as suggested by the simulations of [15], there are 
examples where the It have jumps in time at the limit. So we will only 
be able to prove their right-continuity. 

This regularity in time comes from the stability of the equilibrium defined 
through fll.l2p and Prop. 11.11 Therefore let us define 

/,(t) = /,(/i(Mt, .)})), 

where /j and /i are given by Prop. 11.11 and ip is the uniform limit of cp^ as 
taken in the previous subsection. 
Our first goal is the following result. 

Lemma 3.2 For any fixed s, there exist functions a^, a G C(M+) with 
cr,(0) = a(0) = s.t. 

t 

|/f (r) - Ii{s)\^ dr<{t- s) asit - s) + a{e). 

Remark. Of course the whole point is that and a are uniform in e. It is 
also crucial for the following that a does not depend on s. 

3.3.1 Proof of Lemma [Ml 

Step 0: (f has compact level sets. 

Observe that ips(t = 0,x) — )• — oo when x — )■ ±oo since J^Us(t = 0,x) dx < 
oo and d^^it = 0) is bounded. Because of the uniform convergence of ips{t = 
0) to ip^ on M, one deduces that ip^{x) — )■ — oo when x — )■ ±oo. 
Since d^if e L°°([0,T],M) and hit) E [0, 1], it follows from that dt^p G 
L°°([0, T], M) and thus (p{t, x) -oo when x ±oo for all t > 0. 
Therefore, the set 

Q := {{t, x) G [0, T] X M : ip>{t, x) > -1} 

is compact. 

Step 1: One basic property of {(p = 0}. 

Let us start by the following crucial observation 

Vs, 3ts G C(R+) with Ts{0) = 0, s.t. Vt > s, 

Vx G {(^(t, .) = 0}, 3y G {<^{s, .) = 0} with \y - x\ < Ts{t - s). (3.3) 
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This is a sort of semi-continuity for {(f = 0}. It is proved very simply by 
contradiction. If it were not true, then 

3s, 3ro > 0, 3t„ ^ s, t„ > s, 3yn E {(fitn, •) = 0}, 

diVn, {^{S, .) = 0}) > To, 

where d{y, co) = inf^-g;^ \x — y\ is the usual distance. 

Since all the yn belong to the compact set fl of Step 0, we can extract 
a converging subsequence ?/„ y- As (p is continuous, (p{s, y) = or y E 
{(p{s, .) = 0}. On the other hand one would also have d{y, {(p{s, .) = 0}) > tq 
which is contradictory. 

Step 2: The junctional. 
Denote 

[Is = /i({<^(s, .) = 0}, 
as given by Prop. 11.11 We look at the evolution of 

Fe{t) = I logUe{t,x)dHs{x) = - / i^^{t, x) dns{x) , 

ioY t > s. Compute 



d 
~dt 



Feit) = - I ($^/f(t)r/.(a:)-l ) d[i,{x) + - [ H,{^,{t)) dfi^ 



Now write 



- / Iy2l!it)r]i{x)-l]dij,six) = -^ [ Ue{t,x)dx 

^ JR \i=i J JR 

-- I ( r7i(x) - 1 ) {us{t,x) dx - dfis{x)). 

As E- =1 Vi{^) ~ 1 vanishes on the support of fig, 

- - / I y'(/^(t) - /i(s))r7i(x) ) {ue{t,x)dx -df^six)), 
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with 

m 

Notice that 

k 



A Vi=i / i=i ^i^^>^' 



So we deduce 



[Ueit, x)—Ue{s, x)) dx 




(3.4) 



Step 3: Easy bounds. 
Lemma [3.11 tells that 

-H,{ip,)<CTe. 
The total mass stays bounded in time so 



I {Ue{t,x)—Ue{s,x)) dx < I {Ue{t, x) +Ue{s , x)) dx < C . 

Jr Jr 



And furthermore 



ZL^ (t, x) 1 
log % ' dfis =- I {<^e{t,x) - i^e{s,x))dfls 

Ui;[s,x) e 



If 2 

^ Jr £ 



where the last bound comes from the fact that, by Prop. 11.11 fig is supported 
on {(p{s, .) = 0} C Q, where Q is defined in Step 0. Since in addition we 
know that < 0, 



Jr Ue{S,X) 6 
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Consequently we deduce from fl3.4p the bound 

k 



1 /' ^ Ct^ < e 4 II.. - .|U.,„, . /' ^ (3.5) 



Step 4- Control on A and the measure of {x, ~ 0}. 
For some to be chosen later, decompose 



A{r)dr= j /^(s) ?7j(x) - 1^ Ue{r,x)\^(^r,x)<-a,dxdr 

+ / V'/j(s)?7i(x) - 1 Ue{r,x)I^^(r,x)>-a,dxdr. 
Js Jr J 

For the first part, note again that by fll.9p . there exists R s.t. 

k 

V|X| >i?, < 1/2. 

1=1 

Therefore we may simply dominate 

u^{r, x)I^^(r,x)<-a, dxdr <C (t- s) e~"-'/^ 



Concerning the second part, we constrain 1/2 > > — v'e||L°°(Q) and 
may therefore bound 



,x)>—air 



Ue(r, X)\(r,x)>-2a,- 



Now Yl\=i ^ii^) — 1 is nonpositive on {(p{s, .) = 0} and so 
k \ 

his) r]i - 1 ] I^(^r,.)=o < C sup inf \y - x\ < C r^it - s), 

~[ J a;G{<p(r,.)=0} yG{'P(s.-)=0} 
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by Step 1 as the rji are uniformly Lipschitz. For two sets Oi and O2, define 
in general 

6{Oi,02) = sup inf \x — y\. 

By the same argument, one gets 
^ k \ 

^Ii{s)'r], - 1 j I^(r..)>-2Q, <CTs{t - S) 

+ C6{Mr,.)>-2a,},Mr,.) = 0}). 

Inequality fl3.4p now becomes 

+ C [ T,{r-s)dr + C [ 6{Mr,.)>-2a,},Mr,.) = 0}). 



Conclusion. Eq. (13 ■4p indeed gives Lemma 13.21 if one defines 

1 /•* 

asit - s) = —— / Ts{r-s)ds, 

t S 

a{e) = Ce + 2\\ip, - ^IU-(Q) + CTe-°^/^ 

+ C / 5({y.(r,.)>-2aJ,{¥^(r,.)=0})rfr. 

Of course as is continuous and, as Ts(0) = 0, then trivially o"s(0) = 0. Since 
{^ij, •) > — 2a£} and {v5(r, .) = 0} are subsets of 5"(e) is bounded for 
£ < 1, and thus, in order to complete the proof of Lemma [3.21 we only have 
to check that a{e) — )• when e — )■ for a convenient choice of a^. If we take 
o^e > W^e — Converging to slowly enough to have a^/e +00, we 

only have to prove that 

c! 5({(^(r,.) > -2aJ,{(^(r, .) =0})rfr — ^0 as £ ^ 0. 



By dominated convergence it is enough that for any r 

5({<^(r, .) > -2a J, {<^(r, .) = 0}) 0. 
Just as in Step 1 this is a direct consequence of the continuity of ip. 
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3.4 Compactness of the /f and the obtention of (11.121 ) 

First notice that simply passing to the hmit in Lemma 13.21 



Lemma 3.3 3(t„ G C(M+) with cr,(0) = s.t. Wi 



This means that at any point of Lebesgue continuity of /j, one has li = h. 
We recall that a.e. point is a Lebesgue point for Jj. As the /j were defined 
only almost everywhere anyhow (they are weak-* limits), we may identify /j 
and Jj. This proves f ll.l2p and that Jj is approximately continuous on the 
right for any time t (and not only a.e. t). 

Now let us prove the compactness in L\^^ of each Jf. We apply the usual 
criterion and hence wish to control 







T /-s+h 
h 



Il{t) - Il{s)\dtds. 



Decompose 

"S+h pT 1 ps+h 



j^-j \im-it{s)\dtds< -j^ \im-m\dtds 

j-T i-s+h f-T ps+h 

+ T \k{t)-It{s)\dtds+ - \k{t)-k{s)\dtds. 

Jo h J, Jo h J, 



The first and third terms are bounded directly from Lemmas 13.21 and 13.31 for 
example by Cauchy-Lipschitz 

T ps+h f^/l '"''^^ ^ ^^'^ 



'.yj, _ , 1 , 

''' Js Jo \ 



ii{t)-i,{s)\dtds< I [i- 1 \mt) - ii{s)\^ dt ] ds 

< f {as{h) + d{e)/hf'''ds. 
Jo 
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The second term can be handled the same way after swapping the order of 
integration 

T ^ I'S+h rT+h ^ ft 



It{s) - Ii(t)\dtds = - \It{s) - Ii{t)\dsdt 

Jo Jma.x{0,t-h) 

< / T / \It{s) - ds] dt 

Jo V"' Jmax(0,t-/i) / 

pT+h 

< / {at{h) + a{e)/hy/Ut. 







So finally we bound 

rs+h r'r+h 



- \mt)-ms)\dtds<3 {as{h) + a{e)/hy/^ds 

f-T+h 

<3 / v^(T,(/i) ds + 3(T + h)^a{£)/h. 
Jo 

Since of course the functions (Ts{-) can be chosen uniformly bounded in 
Lemma 13.21 again by dominated convergence, this shows that Vr > 0, 3h, 
3eo{h) s.t. Ve < eo{h) 

/T T rs+h 
- j \It{t)-It{s)\dtds<T. 

This is enough to get compactness of the If in L\^^ and then in any L\^^ with 
j9 < oo, which concludes the proof of Theorem ll.il 
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